A theory of superconductivity in multi-walled carbon nanotubes 
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We devise an approach to describe the electronic instabilities of doped multi-walled nanotubes, 
where each shell has in general a manifold of Fermi points. Our analysis relies on the scale depen- 
dence of the different scattering processes, showing that a pairing instability arises for a large enough 
number of Fermi points as a consequence of their particular geometric arrangement. The instability 
is enhanced by the tunneling of Cooper pairs between nearest shells, giving rise to a transition from 
the Luttinger liquid to a superconducting state in a wide region of the phase diagram. 



Carbon nanotubes offer nowadays a great potential for 
the investigation of mesoscopic physics. Their reduced 
dimensionality leads to very strong electron correlations, 
with the consequent breakdown of the conventional Fermi 
liquid picture[3, H, H, Q. This has become manifest, 
for instance, in the experimental observation of a power- 
law suppression of the conductance over a wide range of 
temperatures j^Q. In a different type of experiments, su- 
perconducting (SC) correlations have been also observed 
in carbon nanotubes at low temperatures. There have 
been measures of a sharp decrease in the resistance of 
ropes with a large number (~ 300) of nanotubes, with 
critical temperatures ~ 1 K 0,0. More recently, quite 
abrupt drops in the resistance have been measured in 
large arrays of multi- walled nanotubes (MWNTs) grown 
in the pores of an alumina template with transition 
temperatures Tc ~ 10 K. 

It has been shown that the growth of the SC corre- 
lations becomes possible in the carbon nanotubes when 
the screened Coulomb interaction is overcome by the ef- 
fective interaction mediated by phonons0,0|. For this 
to happen, it is required in general the coupling between 
a large number of conducting channels, which favors the 
screening of the Coulomb interaction. The case of the 
MWNTs is special in that a certain amount of doping 
seems to be required, to ensure the metallic character of 
all the shells and the possibility of Cooper-pair tunneling 
between them. This is consistent with the fact that the 
drops in the resistance have been only observed in sam- 
ples where all the shells in the MWNTs are contacted by 
the electrodes 0. 

The aim of this paper is to study the low-energy insta- 
bilities of doped MWNTs, which may have a large num- 
ber of subbands crossing the Fermi level in each shell. 
These give rise to a manifold of Fermi points located in 
two disconnected patches around the K and K' points in 
momentum space, with an approximate circular shape 
as represented in Fig. ^ For typical shells with ra- 
dius i? ~ 10 nm, the energy spacing between the sub- 
bands is AE ~ 0.1 eV. Below this energy, the differ- 
ent excitations are given by approximate linear branches 
around each Fermi point, corresponding to the disper- 
sion in the longitudinal momentum p. The electronic 



properties are dictated then by the interaction among a 
collection of one-dimensional (ID) electron liquids, with 
electron fields tpa{p) classified by a Fermi point index a 
that, in our notation, changes sign upon inversion in mo- 
mentum space. In this regime, the scattering processes 
are severely restricted by momentum conservation, and 
the relevant interactions can be read from the hamilto- 
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where there is an implicit sum over all the pairs of Fermi 
points a, b with like chirality (i. c. right- or left-moving 
character). We observe that, as in 2D systems with a 
continuous Fermi line, all the relevant interactions can be 
classified into forward-scattering (/^'[j'' ) , exchange {fj^^b ) ' 

and Cooper-pair icl%^) channels IHIii. 

At the energy scale AE', the main contribution to the 
different couplings in |^ comes from the Coulomb po- 
tential, which in the nanotube geometry is given by 



Vc(r-r') 



with 



[x - x'y + 4i?2 sm^[{y - y')/2R] + 

(2) 

The Fourier trans- 
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form of (O dictates the value of the couplings ^ and 

c'^jj I according to the corresponding momentum-transfer. 
Thus, couplings with momentum-transfer of the order of 
the large momentum Q connecting the K and K' points 
have an intrinsic 1/R dependence, inherent to the be- 
havior of the potential at large momentum. Moreover, in 
the case of the /^"jj^ couplings, the different symmetry of 
the Bloch wave functions for left- and right-moving elec- 
trons leads also to an effective 1/R reduction of their 
strength^, |3, I3- The most important effects of the 
Coulomb repulsion arise in the forward-scattering inter- 
actions, otherwise screened by particle-hole excitations 
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FIG. 1: Schematic representation of the location of the Fermi 
points for the sheUs of a doped MWNT, referred to the hexag- 
onal Brillouin zone (BZ) of graphene. The horizontal lines 
correspond to different low-energy subbands, arising from the 
conical dispersion around the corners of the BZ. The full 
(dashed) hexagon stands for the case of an armchair (zig-zag) 
geometry. The Ca,b {ia,b) couplings label the intratube (inter- 
tube) scattering of Cooper pairs at zero total momentum. 



about each of the Fermi points in the different shells. 
This is consistent with the observation of Luttinger liquid 
behavior and the relative reduction of the exponent for 
the power-law dependence of the conductance in doped 
MWNTs 14, l|. 

At energies much lower than Ai?, however, the cou- 
plings grow large in general, pointing at the onset of 
an electronic instability in the system. The character 
of this instability depends on whether there are differ- 
ent shells electrically active in the MWNTs. Typical 
MWNTs are made of shells with a random mixture of 
helicities. In these conditions, the single-electron tunnel- 
ing between nearest shells is largely suppressed since, in 
general, the misalignment of the respective carbon lat- 
tices leads to a relative rotation and to the mismatch 
of the Fermi points^^, as shown in Fig. ^ However, 
such a mismatch does not prevent the intertube tunnel- 
ing of Cooper pairs with vanishing total momentum. We 
will take this into account by introducing new intershell 
tunneling couplings t^^^ , standing for the amplitude of a 
pair (with electrons around Fermi points a and —a) to 
tunnel into a neighbor shell (with electrons around b and 
—b). To discern between the possible tunneling of singlet 
or triplet pairs, we assign an amplitude t^~^ to the pro- 
cess where the chiralities of the particles in the pair are 



exchanged, different to that for the direct process t'^^^ . 

Below the energy scale Ai5, the electron system is ef- 
fectively ID, and the interaction vertices are corrected by 
diagrams that depend logarithmically on the energy scale 
A measured from the Fermi energy ep- To deal with this 
divergence, one can integrate progressively the electron 
modes starting from the limits of the linear branches at 
£p ± AE Tt'I. This leads to a dependence of the cou- 



plings on the scale I = — log(A/A_B), according to the 
scaling equations 
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where the Vat are defined in terms of the string of Fermi 
velocities Va as Vab — (va + fb)/2. 

The couplings c'^ij are also renormalized by the inter- 
tube tunneling amplitudes and their scaling equa- 
tions turn out to be 
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with h['l^ = 2/^;l, - SaJal 

The above equations apply to the description of 
MWNTs where all the shells participate in the conduc- 
tion, as the t^^ij couplings imply the transport between 
nearest shells. Their relative strength can be estimated 
by taking the square of the interlayer single-electron hop- 
ping in graphite ('--^ 0.1 eV), in units of the hopping rate 
in the graphene lattice (~ 2.5 eV), leading to relative 
amplitudes of the order of 0.002. When the system 
has some tendency to develop superconducting correla- 
tions, such relatively small couplings at the scale AE 
are however reinforced at low energies, together with the 
Cooper-pair couplings. An instability appears then in 
the electron system at scales which are typically in the 
range of Z « 4 — 6. This may be confronted with the 
instance where only one shell is electrically active, as it 
happens in most experiments with MWNTs. Then the 
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^ couplings have to be taken identically zero. The flow 



of Eqs. I©-© still has some unstable regime, but it ap- 
pears at values of the variable Z > 10, which correspond 
to too low temperatures to allow the observation of a new 
phase (given the appearance of Coulomb blockade effects 
in real samples at higher temperatures). 

We have seen that the contribution of the e-e interac- 
tion to the intratube couplings at the scale AE is defined 
exclusively in terms of the potential appropriately 
screened by the RPA sum of finite particle-hole processes 
preserving chirality about each of the Fermi points. To 
complete the evaluation of the couplings, we have added 
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the contribution of the effective phonon-mediated inter- 
action. We have rehed on the results obtained for nan- 
otubes in a wide rang e of diametersfliligt. which coin- 
cide in that the coupling A of the effective interaction for 
small momentum-transfer is about three times smaller (in 
absolute value) than that for momentum about Q [20j. 
We have thus taken a contribution to the original cou- 
plings from phonon-exchange that scales inversely pro- 
portional to the nanotube radius [13 1 running between 
A = -0.1/i? and -0.03/i? (with R expressed in A) for 
the channels with small momentum-transfer, and respec- 
tive values equal to 3A for the channels with momentum- 
transfer about Q 0, ^3 . We have observed anyhow that 
the phase diagram of the system is quite insensitive to the 
variation of A within that range. This is related to the 
fact that, as explained below, the dominant SC instability 
arises essentially from the anisotropy of the interactions, 
created by the peculiar geometry of disconnected patches 
of Fermi points. 

We have characterized the low-energy electronic insta- 
bilities by looking at the growth as A ^ of the different 
response functions, for operators corresponding to order 
parameters for all types of symmetries (charge-density- 
wave (CDW), spin-density- wave, superconductivity with 
various wave symmetries). The SC response functions 
are special in that, for a given (/-wave symmetry, there is 
the intratube response function Rg{uj) and the intertube 
counterpart Tg{uj) measuring the propagation of Cooper 
pairs between nearest shells. The respective derivatives 
with respect to the frequency, Rg = ORg /dloguj and 
Tg = dTn/d\oRUj, satisfy scaling equations in the usual 
fashion 



9 log (A) 

dTgjA) 
a log (A) 



- T.9i'lcii{A)RgiA) 



a.b.s 



E-9i:^a(A)^.(A) 



(8) 
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where ^ is the symmetry factor appropriate for each 
particular wave. We have observed that there is a re- 
gion of the phase diagram, represented in Fig. [3 where 
the growth of the SC response functions is overshadowed 
by the divergence of any of the CDW response functions 
at low energies. Above the phase boundary shown in 
the figure, however, the CDW correlations are small and 
the instability is dictated by the large growth of both 
Rg{uj) and Tg{uj) at low energies. This is the signal that 
the Cooper pairs are able to propagate across the entire 
multi- walled structure. The singularity of Tg{uj) at some 
finite frequency ujc marks the transition to a SC state, 
that we have always found in a channel of p-wave sym- 
metry (together with subdominant d-wave pairing corre- 
lations). 

Quite remarkably, the boundary between the CDW 
and the SC phase in the diagram of Fig. |21 corresponds to 




FIG. 2; Phase diagram of doped MWNTs in terms of the av- 
erage radius R of the shells and the doping level, represented 
by the shift in the Fermi energy ef with respect to the value 
in the undoped system. The upper region corresponds to a 
SC phase with p-wave symmetry. The density plot represents 
in grey scale the values of the logarithm of the transition scale 
uJc, ranging from — log(ti;c/A_E) « 4 (black) to ~ 6 (white). 



the situation where the number of Fermi points changes 
from 12 to 20 in each shell of our model. We have 
solved the scaling equations for a number of up to 50 
Fermi points, which already demands the consideration 
of 1,000 independent couplings. Then, we have found 
that the scale ojc of the SC transition ranges between 
values corresponding to — \og{LUc/ ^E) w 4 — 6. Recall- 
ing that AE ~ 0.1 eV, we see that the predicted critical 
temperatures reach a higher scale of ~ 10 K, which is 
consistent with the experimental observations [9j. 

There is a natural explanation for the appearance of SC 
correlations with p-wave symmetry, which relies on the 
peculiar arrangement of the Fermi points into two circu- 
lar patches. These have a symmetry that corresponds to 
the group C2v ■ The Cooper-pair vertex, taken as a func- 
tion of the angular variables 9 and 9' of the pairs before 
and after the scattering, may be expanded in the form 

V{9, 9') = Vb + Vi cos(6i) cos(6l') + V2 sin(6l) sin(6l') 

cos(26i) cos(26'') + Vi sai{29) sin(26C|0) 

A simple operation allows to estimate the different coef- 
ficients in the expansion. By placing one of the electrons 
in the incoming pair at the maximum angle 9q spanned 
by one circular patch, we may parametrize Vi,V2,V3 
and 14 in terms of the values of the vertex at zero 
momentum-transfer, V{9a,9Q) = J7i, at Q momentum- 
transfer, V{9q^-k ^9q) = U2, with incoming and outgoing 
electrons at opposite angles, V[9q,—9q) = U3, and at 
antipodal points, V{9q,tt + 9o) = U4. We obtain 



Vicos2(0o) 

1^2COs2(0o) 

V4cos2(26'o) 



{Ui -U2 + U3- C/4)/4 (11) 
(f/i + U2-U3- C/4)/4 (12) 
{Ui -U2-U3 + Ua)/A (13) 
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The key point is that, due to the large number of particle- 
hole excitations from the many Fermi points, the inter- 
actions with very small longitudinal momentum-transfer 
are more screened than those with some transverse (as 
opposed to longitudinal) momentum within the same cir- 
cular patch. This means that Ui < U3, so that an attrac- 
tive interaction may exist in the Cooper-pair channels 
with p-wave (V2) and d-wave (V4) symmetry. Moreover, 
as homologous Fermi points in the two circular patches 
can be connected by the fixed vector Q, the screening is 
also enhanced for that particular momentum. This im- 
plies that U2 < U4, leading from Eq. 1^21 to V2 < 0. 
This is the origin of the SC instability, as an attractive 
interaction in one of the pairing channels is always driven 
to strong-coupling in the low-energy limit 13]. 

We stress at this point that the disorder introduced by 
inhomogeneities in the electrostatic potential from the 
different shells or by deformations in a given nanotube 
cannot have a significant effect on the SC correlations. 
A remarkable property of the graphene lattice is that de- 
fects with a spatial size comparable or larger than the 
lattice constant give rise to a negligible backscattering 
in the electron propagation |22j. Then, the effect of such 
inhomogeneities can be modelled by a random potential 
coupled to the uniform density fluctuations ipa^a ■ It 
turns out that the SC correlations are actually insensi- 
tive to this kind of disorder, since the SC order parameter 
only depends on the variable (momentum) conjugate to 
such density fluctuations |2^. In this perspective, only 
lattice defects with spatial size smaller than the lattice 
constant may affect the propagation of the Cooper pairs. 
Such a kind of disorder may prevent or not the develop- 
ment of SC correlations depending on whether it reaches 
some critical strength '2^. This is consistent with the 
experimental observations reported in Ref. 1^, where the 
Raman measurements indicate the absence of defects to 
a large extent in the MWNTs samples. 

We have thus developed an approach that is suitable 
to describe the electronic instabilities of doped MWNTs, 
which may have many subbands at the Fermi level. This 
analysis may be relevant to account for the signatures of 
SC transitions reported in Ref. 0. Our approach relies 
on the scale-dependence of different scattering processes 
at low energies, and it allows to describe the transition 
from the Luttinger liquid to a SC state as some of the 
Cooper-pair couplings start to grow large. 

We have shown that the mechanism of superconductiv- 
ity in doped MWNTs is based on the particular geometry 
of the two patches of Fermi points. A reflection of this 
is the appearance of the SC instability in a channel of 



p-wave symmetry. It is tempting to relate this feature to 
the observation of the interplay between SC and magnetic 
order parameters in some samples of doped graphite 23|- 
All these properties should deserve further investigation 
at the experimental level, for the sake of clarifying as 
well the relation between the superconductivity of differ- 
ent carbon compounds. 
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